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tl 1 Abstract 

In this paper we report the results of a detailed investigation of the complete time evolution 

(N; 

of massless fermion fields propagating in spacetimes of higher dimensional stringy black hole so- 
■ lutions, obtained from intersecting branes in string/M theory. We write the Dirac equation in 

.D-dimensional spacetime in a form suitable to perform a numerical integration of it, and using 
a Prony fitting of the time domain data, we determine the quasinormal frequencies that charac- 
terize the test field evolution at intermediary times. We also present the results obtained for the 
quasinormal frequencies using a sixth order WKB approximation, that are in perfect agreement 

X: 

with the numerical results. The power law exponents that describe the field relaxation at very late 
times are also determined, and we show that they depends upon the dimensionality of space-time, 
and are identical to that associated with the relaxation of boson fields for odd dimensions. The 
dependance of the frequencies and damping factor of the spinor field with the charges of the stringy 
black hole are studied. 
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I. INTRODUCTION 



The study of perturbations around black hole solutions has a long history, beginning 
with the work of Regge and Wheeler in 1957 They addressed the problem of finding 
the effects that the evolution of small metric disturbances can produce on the stability 
of Schwarzschild black holes. To solve the problem, they performed a split of the metric 
perturbation according to the behavior of the different components under the symmetry 
group of the background space-time itself, leading to a classification of disturbances in 
different sectors. The same idea have been used since that time, and a variety of perturbation 
problems has been addressed for a number of interesting solutions, coming from Einstein 
General Theory of Relativity to String Theory 2|-|5]. 

Apart from the obviously important question about the stability of a gravitational system, 
the study of perturbations of black hole solution can serve to know about the geometry of the 
background space-time. It is known that the time evolution of any matter field perturbation 
comes in three stages: an initial transient stage in which the black hole begins to absorb 
the disturbance, and that is highly dependent on the initial form of the fluctuation, followed 
by a set of exponentially damped oscillations called quasinormal modes, and, at very large 
times, late time tails depending of the asymptotic form of the effective potential describing 
the perturbation. 

The above picture is typical for massless perturbations in asymptotically flat space-times, 
where the late time tails have a power law character. For black holes in de Sitter backgrounds, 
the tail depends on the magnitude of the cosmological constant A. For large values of the 
cosmological constant the late time evolution of disturbances is characterized by exponen- 
tial tails, whereas for small A exponential and power law tails can appear. For negatives 
values of cosmological constant, we have usually a time evolution completely dominated by 
quasinormal modes, but in some asymptotically Anti de Sitter (AdS) space-times describing 
charged black holes, there can appear an exponential falloff js S]. 

The quasinormal modes carry information of the black hole space-time geometry, and can 
serve to estimate parameters as the mass, charge and angular momentum of black holes js]. 

Quasinormal modes are also important in such contexts as String/M theory applications, 
in which the celebrated AdS/CFT correspondence can allow us to study the properties of 



strongly interacting systems |lOKL2j|. Related to this general idea, a more specific setup 



2 



called gauge/gravity duality establish that a black hole in asymptotically AdS space-time 
is dual to a strongly interacting system defined at the AdS boundary, and the imaginary 
part of the fundamental quasinormal mode is proportional to the thermalization time in the 
boundary theory [ijj]. Gauge/ gravity correspondence have been applied in to study a 
great variety of physical systems, from quark-gluon plasmas to holographic superconductors 
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16) . Such idea can be extended to systems without Lorentz invariance, as the recently 



discovered asymptotically Lifshitz black holes, giving us the possibility to study strongly 
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18|. 



correlated systems coming from condensed matter theory 

This is the second of a series of papers devoted to the study of perturbations of a family 
of black holes coming from string theory 
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2l| . In this work we extend the results of 



our previous study of four dimensional stringy black holes [22[ to higher dimensional cases. 
We continue considering the evolution of fermion perturbations in the form of massless 
uncharged Dirac fields. 

The motivations to this work coming from the fact that higher dimensional black hole 
solutions appears as natural ones when we consider a String/M theory setup, and for this 
reason it is important to know if the solutions are stable against small perturbations in order 
to make they more reliable. Thus, the propagations of disturbances of different spin weight 
is an important subject to study in such backgrounds. As fermions universally describe some 
matter fields, in most cases they are essential for the structure of the solutions, raising the 



importance to study fermion fluctuations 25|, [26] 



Also the analysis of the quasinormal spectrum of the higher dimensional black holes in 
string/M theory provides a simple way to fix the parameters of the solutions, and conse- 
quently, of string/M Theory itself. Another interesting point is related with the ideas of 
AdS/CFT and holography. In this context, fermion perturbations of the gravitational bulk 
fields equations, turn out to be fundamental matter fields in the boundary conformal field 
theory, providing an important class of objects relevant for condensed matter physics. 

The paper is organized as follows: Section II briefly presents the line element of higher 
dimensional black hole obtained from intersections of branes. In Section III we obtain a 
general wave equation suitable to analyze the propagation of a massless Dirac field in the 
background geometry of the higher dimensional black holes. In Section IV we applied the 
above theory to study numerically the evolution equation of such fields in the considered 
backgrounds, and compute the complex quasinormal frequencies using two methods: Prony 



3 



fitting of time domain data and the sixth order WKB approach. We also presents the results 
of a numerical investigation of the relaxation of Dirac perturbations in higher dimensional 
stringy black holes at very long times. In Section V we consider the large angular momentum 
limit, and write an analytical formula for the quasinormal frequencies. Finally, after the 
concluding remarks, we include an Appendix with Tables showing some numerical results 
obtained in our study. 



II. D-DIMENSIONAL NON-EXTREME STRINGY BLACK HOLE SOLUTIONS 



A particular class of solutions that codify information from string/M-theory and that can 
be regarded as black holes in it, are those non-extremal solutions that appear as a result of 



the intersection of branes in M-theory 
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j|-|2l|. The parameters in those stringy black holes 
solutions can be deduced from the compactification of higher dimensions. The main idea 
of such compactified configurations is to relate static extremal/non-extremal solutions with 
extreme/non-extreme versions of solutions coming from brane surgery. This approach allow 
us to know about the structure of non-extremal black holes, shedding some light about the 



origin of entropy in such systems 



19]. 



is 



In D space-time dimensions, the line that describes the geometry of stringy black holes 



19- 



2l| 



ds 2 = h(ry^f(r)dt 2 + hir^fir^dr 2 + r 2 h(r) 1/D - 2 dQ°- 2 , 



where 



and 



f(r) 



! _ I'ji''" 



h(r) = J] 



i=l 



r J 



(1) 



(2) 



(3) 



The above metrics describes solutions parametrized by r# = [i D - 3 , the position of the event 
horizon, and no charges Qi = sinh 2 (5j), that can be written in terms of the parameters Si, 
coming from the compactification of additional dimensions. The quantity rip depends upon 
the space-time dimensions, being n(4) = 4, n(5) = 3 and n(6 < D < 9) = 2. In the above 
definition of r#, fi is the so called non-extremality parameter, that also contain information 
from the higher compactified dimensions. 
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III. MASSLESS SPINOR PERTURBATION IN SPHERICALLY SYMMETRIC 
SPACE-TIMES 



In this section we write the massless Dirac equation in a general D- dimensional curved 
space-time in a form suitable for its numerical and analytical study. In the following we 
consider a .D-dimensional manifold Ai, whose line element can be written locally as a warped 
product of the form 

M D = O 2 x K°- 2 , (4) 

We introduce coordinates z^ = (y a ,x l ) in Ai D , where y a and x l are coordinates defined in 
the manifolds O 2 and )C D ~ 2 , respectively. In terms of this coordinates, the full metric is 
given by 

ds 2 = Q^dz^dz" = g ab (y a )dy a dy b + w 2 (y a )da 2 D _ 2 , (5) 
where w(y a ) is an arbitrary function defined in O 2 and 

da 2 = r ) i j{x)dx l dx\ (6) 

is the metric on JC D ~ 2 , that we assume is an Einstein manifold. We also consider that J\f m is 
lorentzian. When the line element §5§ describes a black-hole space-time, then the space /C n 
describes the structure of a spatial section of its event horizon. If /C n is a constant curvature 
space with sectional curvature K, then K = represents a flat space, K — 1 a spherical 



manifold and K = — 1 an hyperbolic one 23]. 



In curved space-time the massless Dirac equation is 

= o, (7) 

where y = r M V M is the Dirac operator that acts on the spinor T M are the curved space 
Gamma matrices, and the covariant derivative is defined as = <9 M — ^u^'jalb, with a and 
\i being tangent and space-time indices respectively, related by the basis of orthonormal one 
forms e a = e° . The associated conection one-forms u^ b = u ab obey de a + u a b A e b = 0, and 
the 7° are flat space-time Gamma matrices related with curved-space ones by T M = e M a 7 a . 
They form a Clifford algebra in d dimensions, i e, they satisfy the anti-conmutation relations 
{ 7 «, 7 6 } = -2r] ab , with 7] 00 = -1. 

Under a conformal transformation of the metric of the form : 

= VL 2 g^, (8) 



the spinor ip and the Dirac operator transforms as 241426] 



$ = fH^" 1 ^, (9) 

y^ = iH (D+1) y^. (io) 

If the line element of a space-time takes the form of a sum of independent components 
ds 2 = ds\ + ds 2 ,, where ds\ = g uv (z)dz u dz v and ds\ = g mn {y)dy m dy n then Dirac operator y 
satisfies a direct sum decomposition 

y=y z +y y - (n) 

The above mentioned properties of spinors and the Dirac operator in direct sum metrics and 
conformal related ones, allows a simple treatment of the massless Dirac equations in curved 
backgrounds. For two conformally related metrics, the validity of massless Dirac equation 
in one implies the the validity of the same equation in the other. Then, the idea is to solve 
the Dirac equation in the curved space described by an initial metric tensor performing 
successive conformal transformations that isolate the metric components that depend of a 
given variable, and applied successively the direct sum decomposition of Dirac operator, 
until to obtain an equivalent problem in a space-time of the form M 2 x £ D ~ 2 , where M 2 
is two-dimensional Minkowski space-time and T, D ~ 2 is a (D — 2)-dimensional maximally 
symmetric manifold, in which the s pec trum of massless Dirac operator is known. 



221 ] . First, consider for simplicity a particular class of 



The method work as follows 
metrics of the type (jSJ) 

ds 2 = Q^dz^dz" = -A(y l )(dy ) 2 + B(y l )(dy 1 ) 2 + C{y 1 )dY? D _ 2 {x i ). (12) 

As a basis of one-forms in this space-time we can choose j^] 

e a = e a (y c ) = r b (y c )dy\ a, b =1,2 (13) 
j = e\x k ) = h i j (x k )dx j ) i,j,k = l,...,D-2. (14) 

This basis is convenient because most of the connection one-forms are equal to zero: Uj a b = 
ujj a k = oj a ja = oj a jk = 0. This fact allow us to write for the covariant derivatives V a = 

(2) (2) 

Va <8> I(£)_ 2 ),a = 0,1; V, = Im_?) <g) V a ,j = 2,...,L> - 1, where I(d-2) is the identity 



matrix in (D — 2)-dimensions 27}. Performing a conformal rescaling of the original metric 
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ffT2|) as ds 2 = C{y l )ds 2 , where ds 2 = -£(dy ) 2 + ^(dy 1 ) 2 + dE 2 D _ 2 , then we obtain the 



equation = T^V^V — 0) with ip = C a" ■0. 

In general, we can find a representation of the Gamma matrices in D dimensions 
as a function of their two-dimensional partners: f[ ) 2 ' ) ,f^ 2 ' ) , and the chirality matrix = 

~( 2 ) = p(2)p(2)_ j f D 

is even, the Gamma matrices are of dimensionality 2 a x 2 a , and for 
odd D they have dimension 2~^~ x 2~^~ . We use for odd D the following representation for 
the Gamma matrices 

f { D) = f[, 2) ® I 2 ® I 2 ® I 2 ® • ■ ■ = f[, 2) ® I (c a) 

fJ D) = fS 2) ® I 2 ® I 2 ® I 2 ® • ■ • = ff ) ® I m-aj 

ff } = f (2) g> ® I 2 ® I 2 ® • • • 

ff } = f ( 2 ) ® fS 2) ® I 2 ® I 2 ® • ■ • 

ff> = f( 2 )®f( 2 )®fi 2) ®I 2 ®--- (15) 



fgj = f( 2 )®f( 2 )®f( 2 )®---®fi 2) , 

and for even D we change the last equation above by = f ^ ® T^ 2 ) £g> T^ 2 ^ (g) • • • <g) f g 2 "*. 
It is convenient to define the Gamma matrices T d , d = 2, ...,.D — 1 in the manifold dY? D _ 2 
as: 

r 2 = f f ® i 2 ® i 2 ® i 2 ® • • • 
f 3 = fi 2) ® i 2 ® i 2 ® • • • 
r 4 = f (2) ® fi 2) ® i 2 ® • • • 

f 5 = f (2) ® f[ 2) ® I 2 ® ■ • • (16) 



f J3_i = f (2) ® f (2) g> f (2) (8) • • ■ 
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Then, we can write the last D — 3 curved space Gamma matrices in (fT5|) as: 

tf ) = r i2) ® f 2 

f f > = f ( 2 ) (g) f 3 

ff > = <g) f 4 (17) 



fffi = f< a >®iw 



Also, we have the relations f°( D ) = -T { Q D) , f a ( D ) = fi D) , a = 1, D-l y f ff ( D ) = f { a D) , a = 
2, ...,D — 1. It is interesting to note that under a conformal transformation of the metric 
the Gamma matrices remains equal: = V 1 . 

Using the above representation for the Gamma matrices, the Dirac equation can be 
written as 



(ffoVo + f^Vi) 



I 



(D-2) 



+ r 



(2) 



r l Vi 



-(D-2) 



^ = 0, 



(18) 



where 



is the Dirac operator in the maximally symmetric (D — 2)-dimensional 



space-time, in which we suppose as known its discrete spectrum consisting in the set of 
orthogonal eigen-spinors fa- defined by: 



f fe Vz 



(19) 



being a a set of numbers that label the discrte spectra with eigenvalues k ct = i\ a . Now 
expanding the spinor in the basis 



we can writte the massless Dirac equation in the form 



(20) 



rJ 2) v + rJ 2) Vi±Kf (2) 



# = o, 



(21) 



In the following we work only with the + sign, because the opposite sign case can be worked 
in the same form. We also fix the variables in Af 2 as temporal and radial ones, (t, r). 
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Denoting as ds 2 = — ^pjdt 2 + §^dr 2 the two-dimensional part of ds 2 and performing a 



C(r) 



conformal re-scaling in the form ds 2 = ^ [—dt 2 + dr 2 }, where dr* = \ \dr defines the 



tortoise coordinate we obtain: 



(22) 



being 7* y 7 r * Dirac matrices and 7 s the chirality matrix in the two-dimensional Minkowsky 
space-time with line element ds 2 = —dt 2 + dr 2 . Choosing for the Dirac matrices the repre- 
sentation: 



7* = -ia 3 , 7 r * = a 2 
where the a 1 are the Pauli matrices: 



a 1 = 



1 

1 



a 2 = 



and writing 



, 7 

-i 

1 



^){a 2 ) = -a\ 



a 3 = 



1 
-1 



(23) 



(24) 



we obtain the following equations for the components od the Dirac spinor ip a : 

l ST + &7 + A "*" = ' 



(25) 



where 



dt dr* 



The above equations can be decoupled to obtain 

d 2 (, d 2 („ 



dt 2 



dr 2 



d 2 Xa _ d 2 Xa 
dt 2 dr 2 



+ V + (r)C = 0, 



+ V-(r) X a = 0, 



where: 



V + 



dA <r A 2 

dr* 



(26) 
(27) 

(28) 

(29) 
(30) 

(31) 



This last equations describe the time evolution of Dirac perturbations in curved space- 
times. As the effective potentials V + and V_ are super-symmetric in the sense considered by 
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2 4 6 8 10 



Figure 1: Potential V(r) for Dirac perturbations with I = (bottom) to £ = 5(top). The charges 
parameters of the solution are Q\ = Q2 = Q3 = 1 and the dimension is D = 5 



Chandrasekhar in reference [28], then ( a (t,r) and Xo{t,r) will have similar time evolutions 
and spectra, both quasinormal and scattering. At this point it is interesting to mention 
that for the spinor solution <j!>a , we obtain the same effective potentials. In the following 
we will work with equation (1291 and drop the label (+) for the effective potential, defining 
v} h) (r)=V + (r). 



IV. OBJECT PICTURE, QUASINORMAL MODES AND LATE TIME TAILS 
FOR FERMION PERTURBATIONS 

The general formalism of the previous section can easily be applied to background space- 
times of higher dimensional stringy black holes. The line element ([T]) that describe such 
systems can be easily written as ( !T2|) with the identifications 

A(r) = h(r)~%=*f(r), B{r) = h{r)^ f (r)' 1 , C{r) = r 2 h{r) 1/D - 2 , (32) 

and taking K, n as the (D — 2)-sphere S^ D ~ 2 \ In this specific case in which g?£(£)_2) is the 
line element dQ^ ~ 2 °f the (D — 2)-sphere we have for the eigenvalues of the Dirac equations 



Kg = ±i\e, where = i + ^j^-, £ = 0,1, 2,..., being I the angular multipole number [29]. 



Taking the above elements into account, we obtain for the effective potencial V(r) the 
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t/r H 



10 20 50 100 200 

t/r„ 





Figure 2: Logarithmic plots of the time-domain evolution of £ = (top left), £ = 1 (top right), 
1 = 2 (bottom left) and I = 3 (bottom right) massless Dirac perturbations in (4 + 1) -stringy black 
holes with charges Q\ = 0, Qi = Qs = 1 (top) and Q\ = Q2 = Q3 = 1 (bottom). 



general expression: 



V{r) 



Xtf 



2rh 



(33) 



Figure (CQ) shows the effective potential for different multipole numbers £ in the case of five 
dimensional stringy black holes with Qi = Q2 = Q3 = 1. The form of the effective potential 
is similar for other dimensions and setups for compactification charges. As we can see, the 
effective potential have the form of a definite-positive potential barrier, tending to zero at 
infinity. This character of the effective potential assures the stability of the stringy black 
holes under such spinor perturbations. 

To integrate the equation f j29|) numerically we use the characteristic integration technique 
developed by Gundlach, Price and Pullin [301 . The obtained results can be observed as the 
time-domain profiles showed in Figures (J2]) to (J6j). In such profiles r = 3r# and the time is 
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10 20 30 40 



10 20 50 100 200 

t/r H 



Figure 3: Normal (left) and logarithmic (right) plots of the time-domain evolution of i = 1 massless 
Dirac perturbations in (5 + l)-stringy black holes with charges Q\ = Q2 = 1- 

measured in units of black hole event horizon. 

The first thing that can be observed from this profiles is that the time evolution of 



Dirac perturbations outside D-dimensional stringy b 



ack holes follows the usual dynamics 



for fields in other black hole and brane backgrounds [3l| • After a first transient stage strongly 
dependent on the initial conditions and the point where the wave profile is observed, we see 
the characteristic exponential damping of the perturbations associated with the quasinormal 
ringing, followed by power law tails at asymptotically late times. This picture is independent 
of the chosen values for the compactification charges, and the dimension of the background 
space-time. It its also appreciated that the duration of the quasinormal phase is shorter as 
the space-time dimension is larger. This behavior again is of the same qualitative character 
for different setups of compactification charges. 

To compute the quasinormal frequencies that dominated at intermediate times, we assume 
for the function 0(t,r) in equation (1291) the time dependence: 

Q(t, r) = Z e (r) exp(-iuet) (34) 

Then, the function Z^{r) satisfy the Schrodinger-type equation: 

Zl + [uj 2 - V{r)] Z e (r) = (35) 



d 2 Z t 



drl 

The quasinormal modes are solutions of the equation (I3l)|) with the specific boundary con- 
ditions requiring pure out-going waves at spatial infinity and pure in-coming waves on the 
event horizon. Thus no waves come from infinity or the event horizon. 

12 




Figure 4: Normal (left) and logarithmic (right) plots of the time-domain evolution of i = (top) and 
£ = 3 (bottom) massless Dirac perturbations in (6+l)-stringy black holes with charges Qi = Q2 = 1- 

For the higher dimensional stringy black holes considered in this work, the quasinormal 
frequencies were computed using two different methods. The first method uses directly the 
numerical data obtained previously, and fit this data by superposition of damping exponents. 
This numerical fitting scheme, known as Prony method, allow us to obtain very accurate 
results for the fundamental and sometimes the first overtones . For higher overtones 
is very difficult to be implemented, because we need to do a fitting with a great number of 
exponentials, and also we need to know very well the particular time in which quasinormal 
ringing begins, a difficult point to be solved in general. For this reason we only determined 
the quasinormal frequencies, using this method, for the fundamental overtone. 

The second method that we employed is a semi-analytical approach to solve equation (|35|) 
with the required boundary conditions, based in a WKB-type approximation, that can give 
accurate values of the lowest ( that is longer lived ) quasinormal frequencies, and was used 
in several papers for the determination of quasinormal frequencies in a variety of systems 
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Figure 5: Normal (left) and logarithmic (right) plots of the time-domain evolution of i = (top) and 
£ = 3 (bottom) massless Dirac perturbations in (7+1) -stringy black holes with charges Qi = Q2 = 1- 
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Tables I to X in the Appendix at the end of this report presents the values obtained 
for the quasinormal frequencies with some multipole numbers I for stringy black holes at 
different dimensions and different set of charges. As it is observed, the sixth order WKB 
approach gives results in agreement with those obtained by fitting the numerical integration 
data using Prony technique. 

As it is expected the oscillation frequency increases for higher multipole and fixed overtone 
numbers. Increasing i the magnitude of the negative imaginary part of the fundamental 
overtone (n = 0) increases while for higher overtones the opposite situation arises. For a 
fixed angular number £, the real part of the oscillation frequencies decreases as the overtone 
number increases, and the magnitude of the imaginary part increases. Then, modes with 
higher overtone numbers decay faster than low-lying ones. 

As we expected for stability, all quasinormal frequencies calculated in this work have a 



14 




Figure 6: Normal (left) and logarithmic (right) plots of the time-domain evolution of i = (top) and 
£ = 3 (bottom) massless Dirac perturbations in (8+l)-stringy black holes with charges Qi = Q2 = 1- 

well defined finite negative imaginary part. 

Figure fl7|) shows the dependence of the quasinormal modes with compactification charges 
in the case of five- dimensional stringy black holes. Again the results for higher dimensions 
are exhibit the same qualitative behaviour. As one of the charges increases, the frequency 
and decay factor of the perturbation decreases. For fixed values of the charges, the oscil- 
lation frequency increases as multipole numbers increases, for fixed overtone number. As 
an interesting fact we see that the decay factor of quasinormal modes with n = reaches 
quickly a saturation value for higher multipole numbers. The situation is different for higher 
overtone numbers, as the decay factor of the quasinormal oscillations for fixed n decreases 
as i increases. Again, the qualitative behaviour described is similar for other dimensions. 

Another important point to rise here is the relaxation of the fermion disturbance outside 



the black hole at very late times 42|, |43| . To study the late-time behavior, we numerically 



fit the profile data obtained in the appropriate region of the time domain, to extract the 
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Figure 7: Dependence upon charge Qi, i = 1,2,3 of the real (left) and imaginary (right) parts of 
the massless Dirac quasinormal frecuencies of five dimensional stringly black holes with the other 
two charges fixed Qj = l,j^i. The results correspond to the first overtone for multipolar number 
from I = to I = 4. 

power law exponents that describe the relaxation. 

Figure ([8]) show some representative examples of the results obtained, for different higher 
dimensional stringy black holes and given values of compactification charges. In general, 
all our numerical results indicate that the decay of massless Dirac fields shows a time de- 
pendance in the form of a power law, whose exponent changes not only with the multipole 
number £, but also with the dimension D of the background space-time. Our numerical study 
shows that the power law relaxation factors depends stronlgy of the odd or even nature of 
the space-time dimension D. Stringy black holes of odd dimensions relax at asymptotically 
late times following a power law proportional to t~( 2e+D ~ 2 \ whereas for the case of even 
dimension the falloff is more quickly, following a time evolution proportional to t~^ 2l+iD ^ 
for D = 6, and to r^+sD-ie) for D = 8 

y obtained for the case 



441 ]. In this case it was 



It is important to stress that a similar dependance was analytical 
of scalar perturbations in space-times of Schwarzschild-Tangherlini 
proved that the time evolution at asymptotically late times in odd dimensional space-times is 
independent of the presence of the black hole, because of the Green function of the Laplace- 
Beltrami operator in asymptotically flat backgrounds had always the same tail. For even 
dimensions, the result is strongly dependent on the presence of the black hole and, as in the 
four dimensional case, the tails is determined by the asymptotic behaviour of the effective 
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Figure 8: Late time tails in profiles of time evolution of fermion perturbations in stringy black 
holes. The power-law coefficients were estimated from numerical data represented in the dotted 
line. The full red line is the possible analytical result. The parameters in the figure are D = 5, 
i = 4, Q x = 0, Q 2 = Q 3 = 1 (top left); D = 7, I = 2, Q x = Q 2 = Q 3 = 1 (top right); D = 8, I = 3, 
Q t = Q 2 = Q 3 = 1 (bottom left) and D = 9, 1 = 3, Qi = Q 2 = Q3 = 1. 

potential at larger values of the tortoise coordinate. 



In reference 44[ tails proportional to £~( 2£+3D-8 ) were obtained, for any even dimension 
D. However, as our results shows, fermion perturbations decay following this law only for 
D = 6, whereas the disturbance falloff outside 8-dimensional stringy black holes is more 
slowly. Thus, it seems that the fermion or boson character of the test field have influence in 
the late falloff of the perturbations in even higher dimensions. 

However, we remark at this point that this dependances are only a result consistent with 
our numerical data. A simple analytical argument to support this late time behaviours do 



not exist, in contrast to the case for boson fields, in which the general form o: 



potential is suitable to expand for large values of the tortoise coordinate 44H46| and then 



the effective 



extract the above power law behaviors directly from this asymptotic expansion. The problem 
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related with the analytical determination of the decay factors for fermion perturbations in 
higher dimensional stringy black holes remains open. 



V. LARGE ANGULAR MOMENTUM LIMIT. 

An interesting case arises for large angular multipole numbers because in this limit the 
first order WKB approximation becomes exact and we can obtain an analytical result for the 
quasinormal frequencies. As an interesting fact we should mention that in this limit, for all 
the higher dimensional stringy back holes considered in this report, the effective potentials 
are given by the same general expression 

V(r) = £ 2 A(r) (36) 

where A(r) = J^) r \ ■ Then the first order WKB approximation gives for the quasinormal 
frequencies in this limit the result: 

cu 2 = £ 2 A(r m ) - ii(n + \) >/-2S(r)| r=rm , (37) 

with H(r) = d is given by 

1 



a r 



- {9r 2 f 2 ti 2 + rhf [ti(18f - Wrf ) - Arfh"} 



Ar A fi 

+ 2h 2 [I2f + r 2 f' 2 + 2rf(rf"-5f')}} 

In f )37|) r m the point at which the asymptotic effective potential (136]) reach its peak. This 
value can be determined as the maximum root of the equation 



-f(r) 



2h(r) +r dhi> ' 



+ r h( r )^- = 0. (38) 
dr 



dr- 
ill the particular case of equal charges Qi = Q, the above general equation results 

[(3D - 11)Q + D-1] r 3 - D rg+ 3 - 2{D - 4)Qr 6 - 2D rjP -2r 6 H = 0, D = 5. (39) 

and 

[2(D-4)Q + D-l]r 3 - D r% +3 -(D-5)Qr 6 - 2D r 2 H D -2r% = 0, 6 < D > 9. (40) 
In five dimensions we easily obtain the result: 



r H \Jl + Q + ^1 + Q + Q 2 (41) 
18 



D = 6 (42) 

D = 7 (43) 
D = 8 (44) 

D = 9 (45) 

VI. CONCLUDING REMARKS 

In this paper we studied the evolution of massless Dirac disturbances outside space-times 
corresponding to higher dimensional black hole solutions coming from intersecting branes in 
string theories. 

As in other black hole and brane backgrounds, the time evolution of fermion perturbations 
id divided in three epochs. After the initial transient stage, the evolution is dominated 
by quasinormal oscillations, and at late times by power-low falloffs. We computed the 
quasinormal frequencies for different values of background dimensions and compactification 
charges using two different approaches, a semi- analytical 6th-order WKB formula and Prony 
fitting of the numerical data obtained from the characteristic integration of the evolution 
equations. Both methods gives very close numerical results for the quasinormal frequencies. 
The results for the dependence of the quasinormal frequencies with charges appears to be 
universal for all values of this parameter and all dimensions. 

We also computed the decay factors for the late time relaxation of fermion perturbations, 
showing that for odd space-time dimensions fermion fields relaxes in the same way as scalar 
waves in uncharged backgrounds. However, the late time falloff of the Dirac fluctuatiosn 
for even background dimensions is strongly dependent on the dimension itself, a result that 
contrast with the scalar case. 

Some extensions to this work are interesting. First, it would be an important question 
to determine analytically the late time decay factors, taking into account that for potentials 



whereas for other space-time dimensions we have: 



and 



s/7 /49 O 

r m = r H \l - + 2Q + ^ — + |(8Q + 13), 



6 



r,n =r u \I^Q + 2 + W4+|(25Q + 38), 
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typical of fermion fields do not exist a simple analytical argument to approach this problem, 
as in the case of boson fields. Also, the study of the evolution of boson perturbations: scalar, 
electromagnetic and gravitational ones in this higher dimensional space-times would be an 
interesting subject to be addressed. With the solution of some of those important problems 
we will be involved in future reports. 
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Appendix: Numerical Results for quasinormal frequencies. 



In this Appendix we will show the tables containing the numerical results for the calcu- 
lations of some quasinormal modes of stringy black holes in different dimensions for some 
chosen values of the charges. 



1 


n 


6 th order WKB 


Prony 








0.4349 - 0.1712? 


0.4346 - 0.1709i 


1 





0.7445 - 0.1795? 


0.7443 - 0.1792? 


2 





1.0518 - 0.1806? 


1.0518 - 0.1806? 


2 


1 


1.0073 - 0.5489? 




3 





1.3571 - 0.1809? 


1.3571 - 0.1809 


3 


1 


1.3224 - 0.5469? 




3 


2 


1.2540 - 0.9263? 




4 





1.6623 - 0.1810? 


1.6623-0.1810? 


4 


1 


1.6335 - 0.5458? 




4 


2 


1.5769 - 0.9196? 




4 


3 


1.4942 - 1.3083? 





Table I: Massless Dirac quasinormal frequencies corn for I = to £ = 4. The charges are chosen 
as Ql = 0.5, Q2 = 1, Q3 = 1.8 and the spacetime dimension is D = 5. 
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1 


n 


6* h order WKB 


Prony 








0.5101 - 0.2063? 


0.5098 - 0.2059? 


1 





0.8702 - 0.2152?' 


0.8700 - 0.2151? 


2 





1.2270 - 0.2161? 


1.2270 - 0.2161? 


2 


1 


1.1789 - 0.6579?' 




3 





1.5829 - 0.2163? 


1.5829 - 0.2163? 


3 


1 


1.5445 - 0.6547i 




3 


2 


1.4720 - 1.1106? 




4 





1.9379 - 0.2163?' 


1.9379 - 0.2163i 


4 


1 


1.9063 - 0.6530? 




4 


2 


1.8452 - 1.1014? 




4 


3 


1.7590 - 1.5689? 





Table II: Massless Dirac quasinormal frequencies urn for I = to £ = 4. The charges are chosen 
as Ql = 0, Q2 = 1, Q3 = 1 and the spacetime dimension is D = 5. 



22 



1 


n 


6«« order WKB 


Prony 








0.8545 - 0.3048? 


0.8547 - 0.3048i 


1 





1.2830 - 0.3344? 


1.2831 - 0.3344? 


2 





1.7245 - 0.3365? 


1.7245 - 0.3364? 


2 


1 


1.6107- 1.0275? 




3 





2.1665 - 0.3363? 


2.1665 - 0.3363? 


3 


1 


2.0736 - 1.0214? 




3 


2 


1.8904- 1.7441? 




4 





2.6069 - 0.3364? 


2.6069 - 0.3364? 


4 


1 


2.5292 - 1.0175? 




4 


2 


2.3748 - 1.7245? 




4 


3 


2.1457- 2.4771?' 




5 





3.0462 - 0.3365? 


3.0462 - 0.3365? 


5 


1 


2.9797- 1.0155? 





5 


2 


2.8468 - 1.7130? 




5 


3 


2.6486 - 2.4431? 




5 


4 


2.3878 - 3.2220? 





Table III: Massless Dirac quasinormal frequencies ojrn for £ = to i = 5. The charges are chosen 
as Ql = Q2 = Q3 = 1 and the spacetime dimension is D = 5. 
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1 


n 


6* h order WKB 


Prony 








0.9753 - 0.3387i 


0.9755 - 0.3387i 


1 





1.4500 - 0.3844? 


1.4501 - 0.3844? 


2 





1.9427 - 0.3893? 


1.9427 - 0.3893? 


2 


1 


1.8208 - 1.1924? 




3 





2.4381 - 0.3893? 


2.4381 - 0.3893? 


3 


1 


2.3379 - 1.1858? 




3 


2 


2.1448 - 2.0382? 




4 





2.9323 - 0.3892? 


2.9323 - 0.3892? 


4 


1 


2.8483 - 1.1802? 




4 


2 


2.6839 - 2.0100i 




4 


3 


2.4479 - 2.9088? 




5 





3.4257 - 0.3892? 


3.4257 - 0.3892? 


5 


1 


3.3535 - 1.1768? 





5 


2 


3.2111 - 1.9925? 




5 


3 


3.0037 - 2.8579? 




5 


4 


2.7404 - 3.7965? 





Table IV: Massless Dirac quasinormal frequencies loth for £ = to t = 5. The charges are chosen 
as Ql = 0, Q2 = 1 and the spacetime dimension is D = 6. 
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1 


n 


6* h order WKB 


Prony 








1.2488 - 0.3986? 


1.2490 - 0.3982? 


1 





1.7023 - 0.4570i 


1.7024 - 0.4570i 


2 





2.2068 - 0.4575i 


2.2068 - 0.4575i 


2 


1 


2.0037- 1.4019? 




3 





2.7164 - 0.4547? 


2.7164 - 0.4547? 


3 


1 


2.5460 - 1.3868? 




3 


2 


2.1981 - 2.3835? 




4 





3.2225 - 0.4539? 


3.2225 - 0.4539? 


4 


1 


3.0795 - 1.3759? 




4 


2 


2.7861 - 2.3427? 




4 


3 


2.3271 - 3.3962? 




5 





3.7264 - 0.4539? 


3.7264 - 0.4539? 


5 


1 


3.6035 - 1.3712? 





5 


2 


3.3514-2.3191? 




5 


3 


2.9592 - 3.3271?' 




5 


4 


2.4175 - 4.4380? 





Table V: Massless Dirac quasinormal frequencies urn for i = to i = 5. The charges are chosen 
as Ql = Q2 = 1 and the spacetime dimension is D = 7. 
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1 


n 


6* h order WKB 


Prony 








1.4028 - 0.4138f 


1.4029 - 0.4138?' 


1 





1.8827 - 0.5020? 


1.8829 - 0.5019?' 


2 





2.4272 - 0.5092? 


2.4272 - 0.5092? 


2 


1 


2.2240 - 1.5616? 




3 





2.9827 - 0.5071i 


2.9827 - 0.5071i 


3 


1 


2.8090 - 1.5510? 




3 


2 


2.4626 - 2.6868? 




4 





3.5362 - 0.5061?' 


3.5362 - 0.5061i 


4 


1 


3.3890 - 1.5387? 




4 


2 


3.0927 - 2.6365? 




4 


3 


2.6443 - 3.8594? 




5 





4.0878 - 0.5058? 


4.0878 - 0.5058i 


5 


1 


3.9608,-1.5318? 





5 


2 


3.7039 - 2.6041? 




5 


3 


3.3140 - 3.7652? 




5 


4 


2.7947- 5. 0715i 





Table VI: Massless Dirac quasinormal frequencies loth for i = to t = 5. The charges are chosen 
as Ql = 0, Q2 = 1 and the spacetime dimension is D = 7. 



26 



1 


n 


6«« order WKB 


Prony 








1.6785 - 0.4942? 


1.6786 - 0.4942? 


1 





2.1200 - 0.5864? 


2.1200 - 0.5862? 


2 





2.6741 - 0.5801?' 


2.6741 - 0.5801i 


2 


1 


2.3620 - 1.7865? 




3 





3.2420 - 0.5700? 


3.2420 - 0.5700? 


3 


1 


2.9703 - 1.7521* 




3 


2 


2.4081 - 3.0249? 




4 





3.8019 - 0.5665? 


3.8019 - 0.5665? 


4 


1 


3.5750 - 1.7228? 




4 


2 


3.0967 - 2.9485? 




4 


3 


2.3157 - 4.2988i 




5 





4.3571 - 0.5657i 


4.3571 - 0.5657i 


5 


1 


4.1629 - 1.7103? 





5 


2 


3.7532 - 2.8989? 




5 


3 


3.0894 -4.1785? 




5 


4 


2.1364- 5.6218? 





Table VII: Massless Dirac quasinormal frequencies urn for £ = to i = 5. The charges are chosen 
as Ql = Q2 = 1 and the spacetime dimension is D = 8. 
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1 


n 


6«« order WKB 


Prony 








1.9005 - 0.5183i 


1.9007 - 0.5182i 


1 





2.3570 - 0.6525? 


2.3572 - 0.6525i 


2 





2.9553 - 0.6538? 


2.9553 - 0.6538? 


2 


1 


2.6280 - 2.0151* 




3 





3.5772 - 0.6435? 


3.5772 - 0.6435? 


3 


1 


3.2868 - 1.9856? 




3 


2 


2.6940 - 3.4636? 




4 





4.1926 - 0.6392? 


4.1926 - 0.6392? 


4 


1 


3.9478 - 1.9518? 




4 


2 


3.4387 - 3.3696? 




4 


3 


2.6259 - 4.9856? 




5 





4.8034 - 0.6379i 


4.8034 - 0.6379i 


5 


1 


4.5933 - 1.9350? 





5 


2 


4.1539 - 3.3034? 




5 


3 


3.4546 - 4.8171? 




5 


4 


2.4810 - 6.5801? 





Table VIII: Massless Dirac quasinormal frequencies ujth for I = to I = 5. The charges are 
chosen as Ql = 0, Q2 = 0.75 and the spacetime dimension is D = 9. 
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1 


n 


6«« order WKB 


Prony 








2.1317 - 0.5964i 


2.1317 - 0.596H 


1 





2.5279 - 0.7275i 


2.5279 - 0.7275i 


2 





3.1238 - 0.7074i 


3.1238 - 0.7074i 


2 


1 


2.6867 - 2.1949i 




3 





3.7479 - 0.6832i 


3.7479 - 0.6832i 


3 


1 


3.3458 - 2.1324i 




3 


2 


2.5317- 3.6708i 




4 





4.3556 - 0.6739i 


4.3556 - 0.6739i 


4 


1 


4.0234 - 2.0634i 




4 


2 


3.3140 - 3.5574i 




4 


3 


2.1248 - 5.1437i 




5 





4.9538 - 0.6714i 


4.9538 - 0.6714i 


5 


1 


4.6733 - 2.0331? 





5 


2 


4.0650 - 3.4558i 




5 


3 


3.0514 - 4.9886i 




5 


4 


1.5429 - 6.6723i 





Table IX: Massless Dirac quasinormal frequencies loth for i = to t = 5. The charges are chosen 
as Ql = Q2 = 1 and the spacetime dimension is D = 9. 
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1 


n 


6«« order WKB 


Prony 








2.3363 - 0.5897i 


2.3364 - 0.5897i 


1 





2.7332 - 0.7603i 


2.7333 - 0.7603i 


2 





3.3436 - 0.7571i 


3.3436 - 0.7571i 


2 


1 


2.9325 - 2.3281i 




3 





4.0002 - 0.7351i 


4.0002 - 0.7351i 


3 


1 


3.6023 - 2.2977i 




3 


2 


2.8132 - 3.9808i 




4 





4.6465 - 0.7245i 


4.6465 - 0.7245i 


4 


1 


4.3090 - 2.2305i 




4 


2 


3.6093 - 3.8796i 




4 


3 


2.4595 - 5.6953i 




5 





5.2838 - 0.7208i 


5.2838 - 0.7208i 


5 


1 


4.9972, -2.1939z 





5 


2 


4.3873 - 3.7659i 




5 


3 


3.3968 - 5.5096i 




5 


4 


1.9660 - 7.4862i 





Table X: Massless Dirac quasinormal frequencies ooth for £ = to £ = 5. The charges are chosen 
as Ql = 0, Q2 = 1 and the spacetime dimension is D = 9. 
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